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Abstract 

An expression for the mutual coherence function (MCF) of an electromagnetic beam wave 
propagating through atmospheric turbulence is derived within the confines of the Rytov approximation. It 
is shown that both the first and second Rytov approximations are required. The Rytov MCF is then 
compared to that which issues from the parabolic equation method of strong fluctuation theory. The 
agreement is found to be quite good in the weak fluctuation case. However, an instability is observed for 
the special case of beam wave intensities. The source of the instabilities is identified to be the 
characteristic way beam wave amplitudes are treated within the Rytov method. 

1.0 Introduction 

In the early studies in the late 1950s of electromagnetic wave propagation through atmospheric 
turbulence (Ref. 1 and 2), the weak fluctuation theory known as the Rytov approximation made 
considerable strides in the understanding of the scattering mechanisms inherent in the random permittivity 
field of the troposphere. Within the structure of the theory, propagation quantities such as the log- 
amplitude and phase fluctuations took precedence as well as their characterizing statistical parameters, 
i.e., the corresponding correlation and structure functions of these quantities. Only in the special cases of 
plane and spherical waves were these quantities related to the prevailing statistical parameters of the 
electric field of the wave, i.e., the second order statistical moment known as the mutual coherence 
function (MCF). The case of beam wave propagation was examined (Ref. 3) but an expression for the 
MCF was never given. The theory, however, had its limitations and was replaced with the strong- 
fluctuation theories based on the parabolic equation (Ref. 2). In this treatment, the relevant propagation 
quantities were necessarily the statistical moments of the electric field of the propagating electromagnetic 
wave; the concepts of log-amplitude and phase fluctuations did not, nor needed not, enter into the 
structure of strong fluctuation theory. Furthermore, since the MCF was the statistically most significant 
quantity for applications and experiments, strong fluctuation theory considered not only this quantity for 
the special cases of plane and spherical wave propagation but also for the more general beam wave case. 

It was found, without explanation, that the results for the MCF of the plane and spherical wave cases as 
calculated from the parabolic equation method were identical to those from the Rytov approximation. Due 
to the success of the parabolic equation method, Rytov theory fell by the wayside and the calculation of 
the MCF for the beam wave case was never carried-out within weak-fluctuation theory and a comparison 
was never made to that of strong- fluctuation theory. 

It is the purpose of the present study is to essentially finish bridging the gulf between the Rytov 
approximation and the parabolic equation method by considering the MCF of the beam wave case within 
the Rytov approximation and comparing it to the results of strong fluctuation theory. By definition, the 
weak fluctuation theoretical results should directly follow from those of the strong fluctuation results in 
the limit of weak fluctuations. This will be shown for the beam wave scenario. In the interim, it will also 
be shown that the traditional form of Rytov theory used in the classical studies is incomplete (Ref. 4); the 
subtleties between the order of the Rytov approximation and the order of magnitude of the permittivity 
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fluctuations was overlooked for a complete and unified application of the theory. Finally, a surprising 
limitation of the beam wave modeling within Rytov theory is identified and studied. 

2.0 Recapitulation of the Rytov Approximations 

Using the field decomposition for the electric field of a wave propagating predominantly along the 
x-axis of a coordinate system situated in a random medium, viz., 

E(r) = t/(r)exp(/fcr) (1) 

in the stochastic Helmholtz equation, characterized by the stochastic permittivity field e(r), 

V 2 E(r) + k 2 E(F) + k 2 e(r)E(F) = 0 (2) 


one has for the complex amplitude 

Vlu{f ) + d ^ ) + 2 ik + k 2 e(r)u(r) = 0. (3) 

ox 1 ox 

In the event that the wavelength X of the propagating wave and the size of the smallest inhomogeneity l 0 
characterizing the stochastic permittivity field is such that X « one has that 


d 2 U(r ) 

dx 2 


« 


2 ik 


9 U(r) 


dx 


(4) 


allowing Equation (3) to be approximately written 

Vju(r) + 2ik^^ + k 2 e(r)u(r) = 0 (5) 

ox 

which is a parabolic differential equation of the diffusion type. 

The most straightforward way of dealing with this scenario is to employ Tatarskii’s method of 
solution employing the Rytov transformation (Ref. 2). The idea is to decouple the stochastic factor e(r ) 
from the resulting stochastic field U(r ) in Equation (5). To this end, one can employ the transformation 
(the Rytov transformation) 


U(r) = exp(\|/(r )) 


( 6 ) 


in Equation (5) and obtain 

2?/C ~3x^ + V P ¥ ( f ) + (^ pVt/ ( F ))" + k % r ) = 0 (7) 

which now a non-parametric relation but also a nonlinear one. Solving this equation via a perturbation 
expansion in the quantity v= one has 

t/(r) = U 0 (r)exp(vj/ 1 (r)+vj/ 2 (r)+---) (8) 

in which Uq{? ) is the initial field distribution and 
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¥1 { r ) = -k 2 \G p (r,r%r') -^y d 


is the first Rytov approximation and 


V2(r) = -jG / ,(r,r , )[Vp\|fi(r / )f d 3 "' 


is the second Rytov approximation. Here, the parabolic equation Green function is 


G - {rr)= -\^ 


exp| ik(p— p') 2 /(2(x - x')j 


Writing the complex amplitude U{r ) in terms of an amplitude A(r ) and a phase S(r),i.e., 
U 0 (r) = Ao(r)exp(/5 0 (r)), U(r)= A(r)exp(/S(r)) 

and using these representations in Equation (8) gives 


Vl(r) + V2(r) + - = In -/(So(r)-S(r)) ^ 

= x{r)+iS l (r) 

in which %(r) = ln[A(r)/Ao(r)] is the log-amplitude and Si(r) = S{r) — So(r) is the phase fluctuation. 
Hence, letting ^(r ) = \|/i(r ) + \|/ 2 (r)-l — , one can write for these propagation parameters 

%(r)=Re{\|/ 1 (r)+V|/ 2 (r)+---} = ^( v P(r) + v P !,! (r)) (1 


5i(r)=Im{\|/ 1 (r)+V|/ 2 (r)+---} = ^ 7 ('R(r)- v P !,! (r)) 


These random quantities can only be described by statistical functions such as averages, correlation 
functions, etc. For example, the spatial log-amplitude correlation function is given by 


ByM (Wfl ) + ^(n ))(^2 ) + 'F'fo ))) 

= i (vin yv{r 2 ) + ) + ^(n yv{f 2 ) + )) 

= \ Re{Mfl Mh )) + (^(q )V*(r 2 ))} 

Similarly, for the spatial phase correlation and log-amplitude/phase cross correlation 

B s (n ,f 2 )=\ Re{(^(n Mf 2 )) - (^(n ) x f*(f 2 ))} 
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and 


% (n / 2 ) = \ im{ Wn Mh )) - (^(n )^*(f 2 ))} (i 8) 

However, at this point, it is necessary to consider the relationship between the order in |e| of the statistics 
of V|/, (r) and the order of the Rytov approximation. 

3.0 Connection Between the First and Second Order Rytov 
Approximations— The Necessity of Using Both 

It is assumed, without loss of generality, that the average or statistical mean = 0 . Thus, by 

Equation (9), (vi(F)) = 0. However, by Equation (10), \|/ 2 (^) ~ (e(r))“ and its mean is non-zero. That is, 
one must consider the second order Rytov approximation, which is of second order in the random 
permittivity, for a non-zero mean field. Hence, one has from Equations (14) and (15), to within the second 
order Rytov approximation (i.e., of second order in the permittivity fluctuations) 

(%(F)) = Re{(v|/ 2 (f)}}, (s 1 (F)) = Im{(v|/ 2 (F))} (19) 

Now the correlation functions formed using \|/i (r ) as required by Equations (16) to (18), e.g., 
(\|/i(ri)\|/i(r 2 )) are also of second order in the permittivity fluctuations by Equation (9). Hence, again, to 
within second order in the permittivity fluctuations, only the first order Rytov approximation need be 
used for these quantities. Thus, for a solution to within second order in |e|, by Equations (14) and (15), 
one need only keep the \|/i(F) in the evaluation of the spatial correlation functions. Hence, Equations (14) 
and (15) reduce to 

x(r)=Re{\|/ 1 (r)} = ^(\|/ 1 (r)+V|/t(r)), 5i(r)= Im{^ 1 (r)} = ^-(^(r)- ^i(r)) (20) 

where 

B x {n / 2 ) = ^ Re{(v|/|(f, )\|q (F 2 )) + (\|/i (r, )\\)\(r 2 ))} (21) 

and 

Bs(h ,r 2 ) = ^ Re{(v|/|(r, )\|/i*(F 2 )) - (yi (o Vi ( f 2 ))} (22) 

Similarly, for the cross correlation of log-amplitude and phase fluctuations 

% (h T2 ) = ^ Im{(vi(o Vi (h )} - (vi (n )v* {h ))} (23) 

Therefore, only the product averages (F 2 )\ and (t|/i(fi)\|/i(F 2 )^ need to be evaluated for these 

correlation functions. Of course, from these correlation functions come the associated variances 
o^(r) = By_(f,r), a 2 s (r) = B s (r,r), and a x s(F) = R x5 (f,F). Finally, the structure functions 

D/(a / 2 ) = ((%(n ) - lih )) 2 ) , etc., can be calculated using Equation (20). Both the mean values given by 
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Equation (19) and the correlation functions of Equations (20) to (23) are second order quantities in the 
permittivity fluctuations |e|. 

What is important in the foregoing is the salient fact that the first moment given by the second order 
Rytov approximation, being ~ (e(r ))” , is just as statistically significant as those correlation and structure 
functions found using the first order approximation (Refs. 4 and 5). Any application requiring the 
correlation and structure functions is incomplete without the consideration of the first order moments (%) 

and (Si ) . This circumstance will be demonstrated by the application of Rytov theory for the calculation 
of the MCF of the propagating wave field. 

Complete expressions for the log-amplitude and phase fluctuation statistics of a beam wave are given 
in Appendix A. 


4.0 The Mutual Coherence Function Within the Rytov Approximation 

In the original applications of Rytov theory to stochastic wave propagation, correlation and structure 
functions were calculated and compared with experimental results. The MCF, i.e., the second order 
moment of the associated electric field was relegated to strong fluctuation theory and the parabolic 
equation method. However, with the proper accounting of the second Rytov approximation, an MCF 
expression can be derived based on Rytov theory. 

Using Equation (1), the MCF for a beam wave is defined by 

r 2 (Cpi ,p 2 ) = (E(L, Pi )E*{L,p 2 )) = (u{L, p, )U*(L, p 2 )) ( 24 ) 


Substituting Equation (13) into Equation (8) and using this intermediate result in Equation (24) yields 
r 2 (T,pi ,p 2 ) = Uq(L,Pi )U* 0 {L,p 2 )(exp[x(L,p 1 ) + %(L,p 2 ) + i(S 1 (L,p 1 )- ))] 


(25) 


At this point, it is advantageous to recognize that the ensemble average as indicated in Equation (25) 
can be written in terms of the corresponding characteristic functional 


exp[%(L,pi ) + x{L, p 2 ) + i (Si {L , pi )-S l (L,p l ))] 


= exp 


iq{x{L, pi ) + %(U,p2 ) + i (Si {L,pi ) ■ - Si (L,pj ))} 


(26) 


q——i 


At this point, it is convenient for notational purposes to let %{L,pi ) = %(l), etc., and to write 
X(l) = Xi(l) + X 2 (l) where %i(l) is the first order Rytov approximation for x(l) and % 2 (1) for the second 
order approximation of the quantity; similarly for S(l) = 5) ( I ) + 5' 2 ( 1 ) . Since, as mentioned above, 
quantities up to the second order in the fluctuations |e| are only being considered, one, can perform a 
cumulant expansion of the right side of Equation (26) up to second order in the parameter q and obtain 
(Ref. 6), 


exp 


^{x(^>Pi) + x(i>p2) + '( 5 i( L ’Pi)- 5 i( L ’Pi))}_ 


= exp 


q=-i 


■ exp 


— — q-K 2 H 


K 1 H — K 2 
2 z 


q——i 


(27) 


where the first cumulant K\ is given by 
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K \ - (%i(l) + /C 2 (l) + Xi( 2 )+% 2 ( 2 ) + 7 '{*Si(l)+ S 2 (l)}-/{S 1 ( 2 ) + S 2 ( 2 )}^ 


(28) 


and the second by 

K 2 = ((X 1 (l) + X2(l) + Xl( 2 ) + X2( 2 ) + / ’{^l(l) + 5 2(l)} _; ’{ 5 l( 2 ) + 5 2( 2 )}) ^ 
_ (Xl (l) + X2 (l) + Xl ( 2 )+X2 ( 2 ) + i{Sl (l) + S 2 (l)} - i {^1 ( 2 ) + ^2 ( 2 )}} 


(29) 


From the discussion in Section 3.0, Equation (28) reduces to the simple result 

^i = 2 (X2) (30) 

upon remembering that (% 2 ) and (S 2 ) are independent of transverse position (as indicated by 
Equation (A 13)). Expanding the square terms in Equation (29) and ignoring terms of order higher than 
|e| gives 

2 = ((X(l) + X(2)) 2 )+ 2«((x(l )+ X(2))(5’(l) - S(2))>- ((^(l) - 5(2)) 2 ) (31) 


Finally, substituting Equations (30) and (31) into Equation (27) allows ensemble average in Equation (25) 
to be written 


exp[...]) = exp[i((x(l)+x(2)) 2 ) + i((x(l)+x(2))(yi)-S(2)))-(((s(l)-S(2)) 2 ) + 2(x) 


(32) 


By definition, the phase structure function is 

Ds(l.2>((S(l)-S(2)f) 


(33) 


But 


((xM + x( 2)) 2 } = o|(l)+ 28,(1, 2)+o|(2) (34) 

found by expanding the square within the average. However, again by definition the log-amplitude 
structure function is 


22,(l.2) = ((x(l) - X(2)) 2 ) = cr x(l) - 26,(l,2) + <q(2) (35) 


SO 

((x(l) + x(2)f)=2<T2(l) + 2<r2(2)-D z (l,2) 

Finally, 

((X(l) + x( 2 ))(5(l)-5(2))) = a x ,(l)-%(l,2) + %(2,l)-a )c5 (2) 


(36) 


(37) 
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where it is noted that, in the general case, according to Equations (A9) and (Al 1), B /S (\ ,2) ) 

and c? x s(l) ^ 0^5(2). Using Equations (26), and (27) in Equation (32) gives 


ex P["'])“ ex P 


--^( 1 > 2 ) + Ox( 1 ) + a i( 2 ) + 2 (5C) + / { a %.v( 1 )“ a x.«( 2 ) _ %( 1 ’ 2 ) + %( 2 ’ 1 )} 


( 38 ) 


where D w {\,2) = D x (\,l) + D s (\,l) is the wave structure function. 

An interesting occurrence of a linear combination of the cross log-amplitude/phase correlations 
appears as an imaginary term in Equation (38). Again, for the general beam wave case, this term does not 
vanish. Using Equations (A9) and (All), one has for this linear combination 


G X- ? ( 2 ) %(l>2)+%( 22 ) 

( ul\ L 00 

= -(2 ji ) 2 — J Jlm{ 7 0 (KP( 2 ,l))- 7 o(^(l, 2 ))} exp 


0 0 




O e (k)kc/k<:/x 


(39) 


where P(i , j) = |yp,- - Y*Pj | . This statistical parameter as well as all the other relevant parameters are 

evaluated in Appendix B for the Kolmogorov spectrum of turbulent fluctuations. Substituting these 
results into Equation (38) gives 


exp[- • •]) = exp 


-4 


. 352 £ 2 C,tJ 


L-x 


\5/6 


- 7 / 


Re 1F1 


f 5 kP& V ' 

6” 4(L-x) Y/ 


dx 


7J7 


L (r_ \ 516 

+/(2.176)^ 2 C,7 J 


0 


- 7 / 




hi kF i< ' 

V 6 ” 4(L-x)y /, 


-1^1 


f 5 -p kP n 
V 6 ’’ 4(L-x)y i y 


dx 


(40) 


Equation (40) can be simplified considerably. To this end, one notes that 

P'-(l2)=Pf,=PS 

In addition to this, the integral representation of the confluent hypergeometric function 


1 E 1 (fl,h;z) = 2 — - Ul-tf a '(l + t)" 1 exp(zt/2)<* 
Byci ,b — a ) J 


(41) 


(42) 


yields the fact that 


\F\[a,b-,z*)=\F\(a,b\z ) 


(43) 


for real a and b. Adopting the notation 




kPf 2 


6" 4y,(L-x) 
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one has from the above considerations, 


Thus, it is found that (using Im{z} = (1/2/ )[z - z*] ) 

Im{ i h (P 2 \ )-,T , ( P 2 )} = |[, F{(rf 2 )- lFl ( P 2 ) 

and 

Re{ i Fi (P& )} = |[i^i { p il )+M P & )] 


( 44 ) 


(45) 


(46) 


Substituting Equations (45) and (46) into Equation (40) finally yields the simplified result 


exp[- • •]) = exp 


-4 


.352£ 2 C 2 J[y / (x) 


v V 6 , 

L-x4 „( 5 


in 


kP} 

6’ ’ 4y 7 (x)(L-x) 


dx 


(47) 


As required by Equation (25), one now needs to consider the product of the initial fields of a beam 
wave; using Equation (Al), one has 


WA 


t/ 0 (Z.,pi)t/ 0 - q.P2)=-p-exp| 


<wl 

v iv 2 


{oct (pf + P2 ) ■ + j[oc 2 - (a 2 + al )l]( P 2 - pi )} 


(48) 


where IT 2 = Wq 2 


(l-a 2 T) 2 +cxfL 2 


is the beam radius at a distance L from the output aperture. Finally, 


using equations (47) and (48) in Equation (25) gives for the MCF for a beam wave, to within the second 
Rytov approximation, propagating through Kolmogorov turbulence 


IT 2 

r 2 (L,p 1 ,p 2 ) = ^Lexp| 


•exp 


f W, n 2 ^ 


v w2 ) 


{«i (p 2 + p\ )+ i[ot 2 - (a 2 + a\ )L](p 2 - p \ )} 


-4.352£ 2 C 2 J[y / (x) 


L—x 


^ , 

iFl 


4 ,.; 


kPl 


21 


6” 4y 7 (x)(L-x) 


dx 


(49) 


To make contact with other formulations in what is to follow, it is advantageous to place this result in 
terms of the sum and difference coordinates of pi and p 2 , viz., 

p c = Pl±P2_ 5 p rf = p 2 _pj (50) 

Solving these for p] and p 2 and substituting into Equation (49) and using the definition for P 2 [ gives for 
the beam wave MCF to within the second order Rytov approximation 
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W 2 

r 2 (L,p c ,p rf ) = — exp| 


< wU 

v»U 


2oq 


f n - A 

2 i Pd 


•exp 


-4 


.352£ 2 C 2 jf Y/ (x) 


P? + - 


L- x 


- 2i a 2 - (otf + a? )c (p c • p d 


5/6 


lC 


2 A 

5 %/?Pd-2iYfp c |~ 
6’’ 4 Y/ (x)(L-x) 


( 51 ) 


dx 


5.0 Relationship of the Rytov MCF With That of Strong Fluctuation 
Theory 

It will now be demonstrated that the result of Equation (51) agrees with the weak fluctuation 
approximation of the similar result obtained using the parabolic equation method of strong fluctuation 
theory. First, however, the plane and spherical wave cases will be considered. In the plane wave case one 
has Wo — > 00 and R 0 —> °°. Thus, oq = 0, 0 C 2 = 0, Y i(r|) = 1 and Y 2 (r|) — > 0. In this limit, the asymptotic 
expansion of the confluent hypergeometric function yields 


lim 

Y2^0 


\F\ 


5 p -*P§ 
v 6”4 Y2 (ri)(L-Ti) y 


r(ll/6)l4 Y2 (Ti)(L-Ti) y 


, 5/6 


(52) 


Using this expression in Equation (51) gives for a plane wave 


r2,pH’(Pc ,Pd iL) — r2,pH’(prf vE) — ex P 1 A51 k~C d pJ L 


(53) 


The spherical wave case is defined by W 0 — > 0 and Rq — > Hence, oq — > °°, a 2 = 0, Y i(r|) = r| /L and 

Ya(Tl) — > 0. The same asymptotic form of Equation (52) applies and Equation (51) in this case gives 


1^2 ^ iv ( Pc ,P d )L) ~ 


M 

2 L 


\ 2 


exp 


ik _ _ 
~j~ Pc ■ P d 


exp 


-1.457£ 2 C 2 P ;pJ(q/L) 5/3 


c/q 


f 

V 


kW, 


2\ 


2 L 



ik _ _ 

~j~ Pc ' P d 


ex p[-0.546£ 2 C 2 p^ /3 L] 


(54) 


Considering the case in which p c = 0 (i.e., the MCF about the propagation axis) and normalizing this 
result with respect to r 2iJW (0,0,L) gives 

r 2 jwjiorm (0,p</ ,L) = exp[-0.546/c 2 C|p^ /3 L] (55) 

Thus, Equation (51) gives the same results as does strong fluctuation theory in these limits. 

Now for the comparison to the general beam wave case. Consider the general solution for the beam 
wave MCF as calculated within strong-fluctuation theory, given by Equation (20) to (73) of (Ref. 7), i.e., 

r 2 (p c ,p d ,L) = ^J exp[-np 2 - Z?K 2 + cp d ■ k d + ip c ■ k d - H(p d ,k d ,L)]d 2 K d (56) 


where 
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a = 


2^o 2 


( a!p 

1 + -| 
v a i y 


b= w 2 _ w 0 2 


1 

c = — 

2 


8 8 


(X]L- — (l-a 2 i) 
OC] 


(oqL) 2 + (l-a 2 T) 2 


( 57 ) 


H(p d ,K d ,L) = 2nk 2 J J 


0 — c* 


1 - exps — ik • 


p</ + 




0„(k)(7 2 kJti 


(58) 


Here, O,, (k) is the spatial spectrum of refractive index fluctuations. (The definition of //has been 
slightly modified from that of Ishimaru to facilitate some calculations later in this development.) The 
problem now reduces to showing that Equation (56) reduces to Equation (51) in the limit of weak 
fluctuations. 

5.1 Average Intensity ( p</ = 0) 

At the outset, it is instructive to consider the slightly simpler case of the average intensity of a beam 
wave given in strong fluctuation theory by 


(/(pc £)) = (p c ,0,L) = ^~ J exp[-(?K 2 +/p c • K d - H( 0,K d ,L)\d : 




(59) 


From Equation (58), one has for an isotropic turbulent spectral density 

L < 


H(K d ,L) = 2nk 2 j J 1 - expj —iK • 

0 — - 

L oo 

= (2>')T 2 | j 




V K J 


0„(K)d 2 K</ri 


o o L 


1 -Jc 


IV 

V k j 


<& n (K)KdKdr\ 


(60) 


Using the Kolmogorov spectral density for refractive index fluctuations 

0„(k) = 0.033C 2 k- i1/3 

in Equation (60) and evaluating the resulting the K-integral using analytic continuation yields 

H( K d ,L) = 0.547£ 1/3 C 2 L 8/3 k^ /3 

Since, in the isotropic case, H{yi d ,L) = //(k^,L), Equation (59) becomes in plane polar coordinates 


(61) 


(62) 


(/(p c ,L)) = ^J J exp[-bK 2 !+ip c K d cosQ- H(0,K d ,L)]dQK d dK a 


>2jt 


0 0 


(63) 


W, 


- jexp[-bK 2 - H(K d ,L)]r 0 (K d p c )K d dK a 
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The use of Equation (62) in Equation (63) results in an integral that cannot be analytically evaluated. 
However, in the case of weak fluctuations, i.e., //(k^ ,L) « 1 , the exponential in Equation (63) can be 
expanded to yield 


(l(p c , L) )=^j exp[-bK^](l - H(K d ,L) + ■ ■ -)j 0 (x d p c )K d dK a 


Using Equation (62) in this result allows for an analytic evaluation and gives 




f lib f lV ° (ll o 2 1 

-0.547 r — £ 1/3 C 2 L 8/3 - , F\ — +••• 

U J UJ 6 4b 


Finally, applying the Kummer transformation 


J 11 , p 2 ] f p 2 ) J 5 , p 2 

i F\ — ,1; —exp ] F\ — ,1; — 

6 4 b 4 b { 6 4 b 


and using the definition of b from Equation (57) in Equation (65) results in the average intensity in 
strong-fluctuation theory in the weak- fluctuation limit 

This result should now be compared with that of Equation (51) from Rytov theory for p d = 0; in this 
instance, Equation (51) becomes 

\ V™ '/ W 2 ^ W 2 j 

•exp -4.352^ 2 C,7 J^y 2 (ri)“~^' t^i ~ f >1; dr[ 


Using the definitions of 72 ( 0 ) and ai, the argument of the confluent hypergeometric function gives 


*72(ti)p? ^ 2p| 2p 2 

W 0 2 ot 2 L 2 + (l-oc 2 L) 2 ^ 2 


where the last result comes from using a relationship in Equation (57). Thus, this function becomes 
independent of rj. Equation (68) is then reduced to 
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('(pc^)) = 


*6 

W : 


_ex P 


V w2 J 


•exp 


-4352Fc,? 1 Ff|,l;f|)|(' Y2 (^)^ 


5/6 


dr\ 


(70) 


within which the r| -integration can now be performed upon using, once again, the definitions of 72 ( 0 ) and 
oq as well as Equation (57). This gives 



A-r 1 

k 



j dr \ = 2 516 


"3" 

U, 


£- 5 / 3 ^ 8/3 


/ , \5/3 

V^2 


(71) 


Substituting this into Equation (70) yields the average intensity in the weak- fluctuation limit of Rytov 
theory 


(7(Pc^)) = 


wl 

W 2 


exp 


Ml 

W 2 . 


exp 


-2.9lk ll2 CjjL 812 


( 1 y /3 f 

\F\ 


j 


A i.Ml 

6’ ’W 2 1 


J 


(72) 


the series expansion of which agrees with the result of Equation (67) (this fact was first noted in (Ref. 4)). 
Hence, the average intensity of a beam wave as calculated from strong fluctuation theory agrees with that 
of Rytov theory in the limit of weak fluctuations. 

A digression must now be made with regard to the limits of applicability of the first Rytov 
approximation before graphical depictions of the comparison of general results is given. As derived in 
(Ref. 4), the requirement that the second Rytov approximation can be neglected for correlation functions 
of log-amplitude and phase fluctuations is given by the condition |(%)| « 1 ; for a general beam wave this 
is 


(%} =1.187C ; 7^ 7/6 L 11/6 Re<j(/) 5/6 2 Ti 


5 11 17 iaL ' 

6 6 6 1 + iaL , 


«1 


(73) 


Hence, the limits of applicability are seen to be functions of the parameters of the particular beam wave 
being considered. Take for example, the propagation conditions given by L = 2500 m, C,j = 5x 10 -15 m 2 \ 
and A, = 0.63x10 h m. In the case of a collimated beam wave (R<> — > °°), the behavior of Equation (73) as a 
function of W 0 is displayed in Figure 1. The plot begins at the spherical wave limit and asymptotically 
reaches the plane wave limit. The function |(%)| peaks at the value of Wo ~ v/ZX = .04 m, i.e., that size of 

the first Fresnel zone. As can be seen, the condition of Equation (73) qualitatively breaks down around this 
value. 

Graphical depictions and comparisons of the general results of Equation (51), with p d = 0, and 
Equation (63) are shown in Figures 2 to 5. Here, the normalized intensity of the collimated beam wave 
cases are shown for initial waist radii Wo of 0.0005 m, 0.005 m, 0.05 m, and 0.5 m; the first value is 
essentially for a spherical wave case and the last value effectively represents the plane wave case. In all 
cases, the following propagation parameters prevail: L = 2500 m, A = 0.63 |l, and C 2 = 5 X 10 -15 m- 2/3 . As 
can be seen from Figure 3, the results of Rytov theory begin to diverge at p c = 0.2 m. This phenomenon is 
at its extreme in Figure 4. The results once again coalesce in Figure 5. 
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r 2 (L,o, Pd )/ r 2 (L,o,o) 


Collimated Beam Wave 



0.05 0.1 0.15 0.2 0.25 0.3 

W 0 (meters) 

Figure 1 . — Average log-amplitude for a collimated beam wave. 


Collimated Beam Wave 
W 0 =0.0005 m 



0.1 0.2 0.3 0.4 


p d (meters) 

Figure 2. — Normalized intensity for a collimated beam wave of waist size=0.0005 m. 
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Collimated Beam Wave 
W 0 =0.005 m 


CD 

o 


Q_ 

CD 


(X) = 0.206 


L=2500 m 
?i=0.63 jam 



p d (meters) 

Figure 3. — Normalized intensity for a collimated beam wave of waist size=0.005 m. 


Collimated Beam Wave 
W 0 =0.05 m 


<*>| = 0.521 


L=2500 m 
?i=0.63 pm 



Figure 4. — Normalized intensity for a collimated beam wave of waist size=0.05 m. 
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Collimated Beam Wave 
W 0 =0.5 m 



p d (meters) 

Figure 5. — Normalized intensity for a collimated beam wave of waist size=0.5 m. 


Since the fundamental foundation of Rytov theory and its various approximations are based on the 
parabolic equation, Equation (5), such intensity divergences cannot arise from this fact since they are not 
observed within strong fluctuation theory, which is also based on the parabolic equation. As will be 
shown in Section 5.2, stable results are obtained for the MCF across all ranges of beam wave waist radii. 
In addition, since palatable results for the intensity are obtained for the limiting cases of plane and 
spherical waves, it is thus doubtful that such divergence phenomena are occurring due to the neglect or 
mistreatment of some propagation related mechanism or use of the paraxial approximation. Therefore, the 
beam structure becomes suspect. But again, such is not the case in strong fluctuation theory so whatever 
the source of the intensity divergences are, it is only peculiar to the way the beam wave structure is 
treated within Rytov theory. The study and identification of this intensity divergence behavior will form 
the subject of Section 6.0. 

One can now proceed to examine the complete expression for the MCF as given by Equation (56) and 
compare it to that of the Rytov theory, Equation (51). 

5.2 MCF 

From what has been gleaned from the examination of the intensity, one should consider the weak 
fluctuation limit of Equation (56) from the outset, i.e., 

oo 

W2 . 

r 2 (p c ,p d ,L) = -0- exp \-ap 2 d - bK 2 + cp d ■ K d + i p c • 

071 J 

— OO 

■{\-H(p d ,K d ,L) + ---}d 2 K d (74) 

= if \p c ,p d ^)- r f } (pc >p d,L)+- 
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where 


r 2 ( 1 ) (poPrf^)= -^r ] ex p[~ a Pd~ bK d +c Pd ■ K d+ip c ■ K d ]d 2 K d 


(75) 


and 


rf } (p c ,p d ’ L ) = ^“T J ex Ph«Prf - + c Prf • + ipc ■ Kd]H{pd ,*d ,L)d 2 K d 

— oo 

The evaluation of Equation (75) is straightforward: 

if^Pc ’Prf’ L ) = -^r ex p(-ap</) J exp[/(p c - icp d )-K- d ]exp(-bK2^d 2 K d \ 

— 00 
00 

-exp(-ap 2 ) J / 0 (|p c - icpd \*d )exp(-Z?K 2 )K d dK d 


(76) 


^o 2 


(77) 


(Wo 2 ) 

f 1 ) 

exp 

^ - 
Pc - icpd + 4abp% 

l 4 J 

U*J 

4 b 


Equation (76) becomes, with the use of Equation (58), 

00 00 


r 2 2 ) (poP</^) = -^( 2 ^ 2 ) J J Jexpf-ap^-Z^ + cprf Krf + jpc-Krf] 


—00 0 


I -cxp<! — /K 


P</ + 




V 


(78) 


0„(K)d 2 Krfr|d 2 K</ 


Here, the -integration will be performed first followed by the K-integration. To this end, one has 


W 1 r r 

r 2 2 (p c ,p rf ,L) = ^^ 2 J J{r 1 +r 2 }o„(K)dpd 2 K 


(79) 


—00 0 


where 


T\ = J exp[-up 2 - bK 2 + ?'(p c - icp d )- K d ]d 2 K d 


(80) 


and 


Ti = J exp[-up 2 - 6 k 2 + ?'(p c - icp d ) ■ K rf ]exp 

— 00 

The evaluation of Equation (80) follows that of Equation (77): 


_ _ kti _ 

-IK-Pd-I—Kd 

k 


d 2 


(81) 
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7] = 2n 


f 1 ■) 

exp 

|pc - icp d 

" + 4abp 2 d 

U b) 

4 b 


Similarly, Equation (81) becomes 


Tl = 2n 2 \b exp ( _flp rf _/K ' prf ) exp 


= 271 


Pc — icPd 


Kq 


4 b 


f 1 ) 


\p c -icp d f + 4 abp 2 d 


■ e * fi 2 2 

{2b) 

exp 

4b 

exp 

-iS ■ K r-K- 

4 bk 2 


where 


T1 

-2iy(p c -icpd)+4bp d 


S = 


4 b 


Substituting Equations (82) and (83) into Equation (79) gives 
P 2 (Pc>P</ ’L) = ■ 


K, 

f 1 ) 

\liik 2 exp 

_ Jp£_ 

-icpd ~ + 4abp 2 j 

4 

\2b) 


4 b 


1 1 

0 — 


( 


1-exp 


-/5 k- 


n- „ 2 


V 


4 bk 2 


K- 


0„(K)d 2 K<7ri 


(82) 


(83) 


(84) 


(85) 


Using the isotropic Kolmogorov spectral density given by Equation (61) yields for the K-integration 


( 


1-exp 


-/5 k- 


,.2 


V 


4 bk 2 


K- 


;O i; (K)d 2 K = 0.033(27t)C ; 7 

K _8/3 <7k 


J '-#) 

0 - 

= (2tx)(0.033)C i 


exp 


p —2 


4 bk 2 


K- 


( 86 ) 



( 

— F 

f--T 

T q 2 ^ 

5/6 

, F, 

f 


,_i2 A 

5 k 2 b 

,2, 

i 

V 

V 6 >y 

[4bk 2 ) 


V 

6 ’ 

’ fi 2 

) 


where analytic continuation was used to obtain the final result. Using this in Equation (85) and 
substituting Equations (77) and (85) into Equation (74) finally yields 
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( 87 ) 


r 2(Pc»P</>^) 



f 1 ') 

exp 

Pc 

1 

C4V3 

CL 

^t 

+ 

<N 

•73 

1CL 

.O 

T 

4 

\ 2 b) 

4 b 


r 


L 



V 


0 


iV 


\5/6 


4M 2 


|V, 



sf k 2 b 

r\ 2 


\ \ 

dr\ + • • • 

7 


Using the definitions for a, b, c, and S, and employing the transcription r) -» L-r\, one obtains, after 
some rather involved algebra, 


(pc >P cl >V) — 


(a x L) +(l-a 2 )‘ 


exp 


ai 


< 0 2 4 

p?+^ 


■ ?(a 2 - (oif + oc 2 )z,)p c • 


(a 1 L)“ + (l-a 2 ) 


( 88 ) 


1- 


435 \k 2 C 2 ) f-y 2 (ri) 


, V-fi 


\5/6 


1G 


= A 


_5 fc|Yiprf~ 2 tY 2 p c 
6” 4y 2 (L-ri) 


c/q + • 


The first two factors of this expression describe the overall beam wave structure; what is important to 
note here, however, is that the expansion given by the third term agrees with the corresponding weak 
fluctuation expansion of Equation (51). 

Graphs of the normalized MCF’s of the entire expression of Equation (51) with that of Equation (56) 
are shown in Figures 6 to 8 for a collimated beam wave subtending the spherical to plane wave limits. 
The largest discrepancy between these results occurs at W 0 = 0.05 m; however, the Rytov results do not 
diverge as they do for the associated intensity of Figure 4. 


Collimated Beam Wave 
W 0 =0.005 m 



p d (meters) 

Figure 6. — Normalized MCF for a collimated beam wave of waist size=0.005 m. 
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Collimated Beam Wave 



0.005 0.01 0.015 0.02 0.025 0.03 


p d (meters) 

Figure 8. — Normalized MCF for a collimated beam wave of waist size=0.5 m. 
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6.0 Limits of Beam Wave Parameters Imposed by Their Use in the Rytov 
Approximation 

As shown in the last section, the calculation within the first and second Rytov approximations of the 
intensity of a beam wave propagating through turbulence revealed a divergent instability in the 
predictions for beam waves of waist radius on the order of the prevailing Fresnel zone. The MCF, 
however, remained stable in this situation. It thus becomes important to identify the source of this 
divergence for such beam waves. 

As shown in previous treatments, all orders of approximation of Rytov theory possess the 
combination of initial field distributions 


where for a Gaussian beam wave, one has Equations (Al) and (A2), viz. 


U 0 {x,p') = 


_L 

1 + z'ax 


exp 


( ka\ p' 1 


V * J 


1 + iCLX 


with 


2 

OC = OC ] + ZOC? » OC ] = — 

kWg 



(89) 


(90) 


(91) 


For simplicity, and without loss of generality, consider the collimated beam wave case in which Rq — > <=°, 
i.e., a.2 =0. Substituting Equation (90) into Equation (89) gives 


where 


R(L, p;x,p') 


1 + i ct i L 
1 + z'OCiX 


exp 


ka.\ | p' 2 

2 [l + z'ajx 


P 2 

1 + i(X\L 


1 + i (X i L 
l + /aix 


exp 




1 + /oqL 
l + /aix 


exp 


k 

2 




1 + z ex i L 

exp 

1 + z'ajx 


ik 

A+-S 

2 



(92) 


( 93 ) 


is the amplitude factor describing the evolution of the wave amplitude and 
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X 


L 

\2 


(94) 


S = 


V2. 


a i 


+ x- 


a. , 


+ Lr 


is the phase factor. 

An important condition for the convergence of any integral expression in which Equation (92) finds 
itself is A < 0. Thus, consider the amplitude term within the exponential of Equation (92) more closely, 
i.e., 


a^-\ K ( x )p' 2 +\k(l) p 2 , 



(95) 


Look now at the first term of on the right side of Equation (95) and consider the two cases (l/ai )” » x 2 
and ( 1 / a | )“ « x 2 , i.e., using the definition of ai, Wq > -Jkx and Wq < V^x . 


6.1 Case (i) (l/ai )“ » x 2 

Take the worst case in which x = L. Thus, A'(x) ~ K(l) = aj and one has for the first term on 
Equation (95) 

~^K(x)p' 2 ~ ^K(L)p' 2 = ~ a ip /2 = ~j (96) 

For any integral over p' in which the expression of Equation (92) occurs, most contributions over the 
integration range come from p' 2 /u / ( j 2 < 1 . Taking the equality to hold in the worst case and using this 
condition in Equation (96) yields 

^(x)p' 2 ~^(L)p' 2 ~l (97) 

Substituting this intermediate result into Equation (95) gives for the condition of convergence of any 
integral involving Equation (92), i.e., A < 0, 

A = -l + ^K(L)p 2 <0 (98) 


6.2 Case (ii) (l/oci)“ « x 2 


Once again, taking x = L, one has K[x) = (l/cii )/L 2 . Hence, for the first term on the right of 
Equation (95) 

k 1/a! „, 2 _ k 2 W 2 




■P - 


4L 2 


(Wo 2 ) 

f p,2 l 

(uj 

[xlJ 


(99) 
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But by assumption, W 0 2 « A L . However, in most propagation scenarios, p' > VAL in the significant 
range of integration over p'; there is thus an indeterminate size for the last term of Equation (99). At this 
point, without further analysis, it will be taken to be ~1. Therefore, once again the condition 


f*A)l p'2-|k(L)p' 2 -1 
obtains and Equation (95) gives for a convergent integral condition 

A = -l + ^K(L)p 2 <0 

which is identical to the result in Equation (98) of Case (i) (Sec. 6.1). 
Using the definitions given above, Equations (98) and (101) reduce to 


rax | 

1 + afZr 


<1 


which can be further reduced to 


P 2 /w 2 

l + A 2 L 2 /7t 2 W 0 4 


This finally can be simplified once more to yield the constraint 


p 2 < W 0 2 + 


y:-c- 

71 2 Wq 2 


( 100 ) 


( 101 ) 


( 102 ) 


(103) 


(104) 


existing between the beam wave and propagation quantities. 

This condition holds for all types of collimated beam waves and can also be a good indicator for 
convergent and divergent beam waves. It gives the possible range of values of the transverse coordinate 
for which the intensity of a beam wave remains stable within Rytov theory. For plane (Wo — > °°) and 
spherical waves (W 0 — > 0), the condition holds trivially. The condition shows itself for nominal values of 
Wo. For example, as shown in Figure 3, a case in which the intensity divergence was noted to occur is 
given by W 0 = 0.005 m, A =0.63 pm and L = 2500 m. This gives through Equation (104) p < 0.1 m, which 
is just before the intensity divergence sets in. However, the case for which W 0 = 0.05 m shown in 
Figure 4, i.e., right at the Fresnel zone length, gives p < 0.05 m; the intensity divergence occurred well 
before this value. The discrepancy is be due to the fact that an asymptotic analysis on either side of the 
quantity VA L was used to derive Equation (104). One can expect a disagreement in situations in which 
Wo ~ VAL. A more complete analysis than that given above is needed to cover such cases. Suffice it to 
say that the condition Wo ~ VA L precludes the use of Rytov theory to describe the behavior of beam 
wave intensity. 

Thus the source of intensity divergences for beam waves within the confines of Rytov theory has 
been identified and an expression is given for the stable range of intensity predictions. The fact that the 
MCF remains stable over these transverse ranges is connected with the fact that the difference coordinate 
over which the MCF is evaluated is symmetric about the beam wave axis; the centroid coordinate is taken 
to be zero in these cases. Thus, the divergences that arise off axis are more or less canceled by performing 
the field comparison across the difference coordinates as is done to obtain the MCF. 
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7.0 Conclusion 


A closed form solution for the MCF based on the first and second Rytov approximations has been 
derived and compared with that from strong fluctuation theory. The agreement between these two 
approaches is quite good in all beam wave cases that satisfy the weak fluctuation case save for the 
intensity distribution for a collimated beam wave. In the region where the initial beam waist size is on the 
order of the first Fresnel zone length, the solution in the Rytov case diverges. The source of this 
divergence is found to be in the way a beam wave is modeled within the Rytov approach; it arises from 
the ratio of the complex amplitudes that occurs within the theory. 
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Appendix A 

One has for a beam wave of waist radius Wq and a phase front radius of curvature R 0 both measured 
at the output aperture of the transmitter, an initial field distribution of unit amplitude given by (Refs. 2 
and 3) 


Uo{r) = U Q {x, p) 


1 

exp 

1 + iax 


' ka N 

l~2~. 


P 2 

1 + iax 


where 


2 

(X = OCi + iOC? ? OC ] = 7T 

kwj 



using Equations (9) and (11) with Equation (Al), one obtains the correlation expressions 


(Vi(APi)Vi (AP2 



o 


J exp(/K-g)// 2 (L- a',k)F £ (k)J 2 kJx 


(Al) 


(A2) 


(A3) 


and 


t|/ 1 (L,p 1 )\|/t(L,p 2 )) = ^-J jexp(iK-P)\H(L-x,K)\ 2 F E (K)d 2 Kdx 


0 — 


where the complex coordinates are Q = y(pi - p 2 ) and P = ypj — y ¥ p 2 with 


/ T \ 1 + iax . 

y = y (L,x ) = —— — = y R ~iyi, Y/>0 

1 + iaL 

(l-a 2 x)(l-a 2 L) + afLx _ a^L-x) 
(l-a 2 L) 2 +afL 2 ’ (l-a 2 L) 2 + afL 2 


(A4) 


(A5) 


and 


H(L- x,k) = exp 


7K- 

~2k 


(L-x)y 


(A6) 


Finally, the two dimensional spectral amplitude of the permittivity fluctuations is given by 
F e ( k) = 27iO £ (k) where 0 £ (k) is the associated three-dimensional spectral density. 

The solution to the second Rytov approximation yields (Ref. 4) 

Lx 00 

<|i|/ 2 (L,p)^ = ^\|/ 2 (l)^ = -/ — J J J y 2 (x,x')// 2 (x- x / ,K)K 2 F e (K)d 2 KJx / dx (A7) 

0 0 — 

a result that is independent of the transverse position. 

From Equations (21) to (23), one obtains in the case of an isotropic spatial spectrum of fluctuations 
(Ref. 3) 
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and 


In particular, 


and 


B x (L,pi,p 2 ) = (2jt) 2 ^-IJ J Re<j/ 0 (i^)exp 


o o 


-\{ L ~ x )v 


+./ 0 (xQ)cxp 




>O f (lc)Kc/Kc/x 


Pl>p 2 ) = -(2tt) 2 J Im<j/o(icP)exp 


v L OO 


+/ 0 ( K Q)exp 


o o 


~( L -^)y 


-^-a>y / 


dX.fKjKc/Kc/x 


( k 2 ' 1 


0 0 


a 2 (L,p) = (2jt) 2 — |J J Re|/ o (2iy. / Kp) exp 

lo e ( K ) KclKdx 


K 


Y( L ~ x h 


+ exp 


-\{L-xb 


hsi L ’ p) = -(2jt)^ 


( k 2 


|j J Imj J o (2 i'Y/Pk) exp 


o o 


~\{L-x)y I 


+exp 




(l> f: (K)KdKdx 


In addition, one has for the log-amplitude and phase structure functions 

D % (L,pi,p 2 ) = (2jt) 2 J Re { ( J o ( 2? Y/ k Pi ) + J() (2iy /Kp 2 ) - 2J 0 ( kR)) 


0 0 


•exp 




+ 


2(l-7 0 (Kg))exp 




>O e (K)K(:/Kc/A' 


Finally, for the mean log-amplitude and phase fluctuations (Ref. 4) 


Md)L 

h(i))l 


7 2K) Um [?)I I I A«')exp 


L x 


0 0 0 


x')y(x,x') 


k 2 0 £ (k)k dKdx'dx 


(A8) 


(A9) 


(A 10) 


(All) 


(A 12) 


(A 13) 


In these expressions, y / = — Im{y} >0. Before an evaluation of these various equations is given for the 
Kolmogorov spectrum of turbulence, an expression for the MCF will be developed. 
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Appendix B 

Using the Kolmogorov spectrum for permittivity fluctuations O e (k) = 0.033Cf K -I1/3 where, in terms 

of the refractive index structure constant C| = 4 C 2 , and evaluating the requisite functions needed in 
Equation (38), one obtains the following results from Equations (39), (AlO), (A12), and (A13) 


Av(l>2) = D W /(L,p 1 ,p 2 )- D x (L,p! ,p 2 )+D s (L,pi ,p 2 ) 


= -0.816C 2 £ 7/6 L 11/6 


v 5/6 


a X L 


( f 
iFi ~ 

V v 


5 , 2p( 

6’ ’W 2 


2 ) 


(l-a 2 L) +a^L 2 
■>2 h 


Ml 


Rff +4.352*^1 

'J 0 


L — x 


\ 5/6 


Y / 


•Re lF, 


' 5 *P,j j V 

v 6 ” 4 ( l -^)y/, 


r/v 


a 2 (L,p) = 2476R /6 L 11/6 C 2 


r r 

( 6 ) 

( 

Rd 

OR 

\i ‘- 5/6 2^1 

V 


5 11 17 iaL 


\ 


6 6 ’6 1 + iaL 

2 h 


"3" 

J, 


,5/6 


aiL 


(l-a 2 L) 2 -t-a^L 2 
(x) = -2.176fc 7/6 C 2 L I1/6 Re 


lU 


_5 2p 

, 4? 


6 ’ ’ IE 2 


v l ly 


^\ 5 /6 { 5 11 17. /«L 

V 6 6 6 ’ 1 + iaL j 


A %5 =2.176£ 2 C 2 J 


N 

1R 

i 

• — 

4 

5/6 \ f 

IiJ iFi 

0 V ~ J 

{ v 


4 , 1 ;- 


kP^ x 


\ 


4(L-a-)y/ 


-iFj 


' J_ tf , 2 2 A 

V 6" 4(Z.-j)y, 




(Bl) 


(B2) 


(B3) 


(B4) 


NAS A/TM— 20 11-216974 


27 



References 


1. V.I. Tatarskii, Wave Propagation in a Turbulent Medium, McGraw-Hill, New York, 1961. 

2. V.I. Tatarskii, The Effects of the Turbulent Atmosphere on Wave Propagation, U.S. Dept, of 
Commerce, TT-68-50464, Springfield, VA, 1 97 1 . 

3. A. Ishimaru, “Fluctuations of a Beam Wave Propagating Through a Locally Homogeneous Medium,” 
Radio Sci. 4 (4), pp. 295-305 (1969). 

4. R.M. Manning, “Beam Wave Propagation within the Second Rytov Perturbation Approximation,” 
Radiophysics and Quantum Electron. 39 (4), pp. 287-295 (1996). 

5. R.M. Manning, “Theoretical Investigation of Millimeter Wave Propagation Through a Clear 
Atmosphere,” Proc. SPIE 410 , Laser Beam Propagation in the Atmosphere, J.C. Leader, ed., 
pp. 119-136 (1983). 

6. A.D. Wheelon, “Skewed Distribution of Irradiance Predicted by the Second-Order Rytov 
Approximation,” JOSA A 18 (1 1), pp. 2789-2798 (2001). 

7. A. Ishimaru, Wave Propagation and Scattering in Random Media, Academic Press, New York, NY, 
Vol. 2, 1978. 


NAS A/TM— 20 11-216974 


28 



REPORT DOCUMENTATION PAGE 

Form Approved 
OMB No. 0704-0188 

The public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, gathering and maintaining the 
data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this collection of information, including suggestions for reducing this 
burden, to Department of Defense, Washington Headquarters Services, Directorate for Information Operations and Reports (0704-0188), 1215 Jefferson Davis Highway, Suite 1204, Arlington, VA 22202-4302. 
Respondents should be aware that notwithstanding any other provision of law, no person shall be subject to any penalty for failing to comply with a collection of information if it does not display a currently valid OMB 
control number. 

PLEASE DO NOT RETURN YOUR FORM TO THE ABOVE ADDRESS. 

1. REPORT DATE (DD-MM-YYYY) 
01-01-2011 

2. REPORT TYPE 

Technical Memorandum 

3. DATES COVERED (From - To) 

4. TITLE AND SUBTITLE 

The Atmospheric Mutual Coherence Function From the First and Second Rytov 
Approximations and Its Comparison to That of Strong Fluctuation Theory 

5a. CONTRACT NUMBER 

5b. GRANT NUMBER 



5c. PROGRAM ELEMENT NUMBER 

6. AUTHOR(S) 

Manning, Robert, M. 

5d. PROJECT NUMBER 


5e. TASK NUMBER 


5f. WORK UNIT NUMBER 

WBS 182306.01.02.03 

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
John H. Glenn Research Center at Lewis Field 
Cleveland, Ohio 44135-3191 

8. PERFORMING ORGANIZATION 
REPORT NUMBER 

E- 17615 

9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 

10. SPONSORING/MONITOR'S 
ACRONYM(S) 

NASA 


11. SPONSORING/MONITORING 
REPORT NUMBER 

NAS A/TM-20 11-216974 

12. DISTRIBUTION/AVAILABILITY STATEMENT 

Unclassified-Unlimited 

Subject Categories: 32, 43, and 70 

Available electronically at http://gltrs.grc.nasa.gov 

This publication is available from the NASA Center for AeroSpace Information, 443-757-5802 

13. SUPPLEMENTARY NOTES 


14. ABSTRACT 


An expression for the mutual coherence function (MCF) of an electromagnetic beam wave propagating through atmospheric turbulence is 
derived within the confines of the Rytov approximation. It is shown that both the first and second Rytov approximations are required. The 
Rytov MCF is then compared to that which issues from the parabolic equation method of strong fluctuation theory. The agreement is found 
to be quite good in the weak fluctuation case. However, an instability is observed for the special case of beam wave intensities. The source 

of the instabilities is identified to be the characteristic way beam wave amplitudes are treated within the Rytov method. 

15. SUBJECT TERMS 

Electromagnetic wave transmission; Wave propagation; Coherent scattering; Coherent radiation 


16. SECURITY CLASSIFICATION OF: 


17. LIMITATION OF 

18. NUMBER 

19a. NAME OF RESPONSIBLE PERSON 




ABSTRACT 

OF 

ST1 Help Desk (email:help@sti. nasa.gov) 

a. REPORT 

u 

b. ABSTRACT 

u 

c. THIS 
PAGE 

u 

uu 

PAGES 

34 

19b. TELEPHONE NUMBER (include area code) 

443-757-5802 


Standard Form 298 (Rev. 8-98) 
Prescribed by ANSI Std. Z39-18 




































